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We demonstrate a five-bit nuclear-magnetic-resonance 
quantum computer that distinguishes among various func- 
tions on four bits, making use of quantum parallelism. Its 
construction draws on the recognition of the sufficiency of 
linear coupling along a chain of nuclear spins, the synthesis of 
a suitably coupled molecule, and the use of a multi-channel 
spectrometer. 



I. INTRODUCTION 

While quantum computers of two bits have been imple- 
mented jy, as have nuclear-magnetic-resonance (NMR) 
quantum computers of three bits Q| , extending the num- 
ber of bits has not proved easy. We report the implemen- 
tation of an NMR quantum computer having five bits, in- 
volving the use of a linear coupling pattern (3) , synthesis 
of a molecule having five usable spin-active nuclei with 
predominantly linear spin-spin coupling, and the devel- 
opment of radio-frequency (r.f.) pulse sequences to act as 
quantum logic gates for the molecule synthesized. Tech- 
niques to suppress unwanted couplings between nuclear 
spins are described, as are techniques to avoid perturb- 
ing some nuclear spins while manipulating others. Re- 
sults are presented of a test of the five-bit computer on 
a problem of Deutsch and Jozsa to distinguish one class 
of mathematical function from another W. 



II. DEFINITION OF AN n-BIT NMR COMPUTER 

An ro-bit quantum computer is called on to do three 
things: 1) accept an instruction to prepare a starting 



state and prepare that state; 2) accept instructions for 
and implement quantum gates (from which more general 
unitary transformations of the state can be composed); 
and 3) measure the state and yield an outcome. The 
connection to computation with classical computers de- 
pends on the recognition, due to Bennett ||, that all 
classical computations can be made reversible. Any ter- 
minating reversible computation is a permutation of the 
inputs, which is unitary, and thus belongs to the class 
of transformation performable on a quantum computer. 
(For issues of possibly nonterminating programs, see ||.) 

In theory, a variant of the quantum computer is the 
expectation- value quantum computer (EVQC), which in 
place of an outcome of a measurement yields the expec- 
tation value |t],|8| . NMR quantum computing was born 
of the recognition that an EVQC can be approximated 
by use of an NMR spectrometer containing a liquid sam- 
ple, the molecules of which have n atoms with a nuclear 
spin of 1/2 (and possibly other atoms, either spinless or 
having spins not used) ■ Because tumbling of the 

molecules decouples each molecule from all the others, 
the sample can be 

described by a density matrix for the nuclear spins of 
the atoms of a single molecule |l(J, with only the spin- 
degrees of freedom, corresponding to the desired Hilbert 
space of dimension 2™. NMR spectrometers sense only 
the traceless part of the density matrix, so in place of 
matter in a pure state, an NMR computer can use a 
liquid sample described by a density matrix proportional 
to a sum of a pure state and any multiple of the unit 
matrix. Such a density matrix, called a pseudopure state 
, plays a role in the 5-bit quantum computer. 

Acting as an n-bit EVQC, a suitable NMR spectrome- 
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ter allows the preparation of a pseudopure starting state, 
the programming and execution of r.f. pulse sequences 
that implement quantum gates, and the 

determination of expectation values visible in NMR 
spectra. To perform the unitary operations required of 
a quantum computer, a sufficient set of quantum gates 
consists of all single-spin operations and all controlled- 
not gates that act on one nuclear spin under the con- 
trol of another nuclear spin. Single-spin gates are imple- 
mented by selective r.f. pulses. Controlled-not gates be- 
tween nuclei having spin-spin coupling will be described, 
along with techniques to avoid unwanted influences on 
other spins. A key feature of the present design of the 
NMR quantum computer is the reliance on a chain of 
linear coupling and the use of swap gates to implement a 
controlled-not in which a spin j controls spin k, where j 
and k have no direct spin-spin coupling H. This allows 
use in NMR quantum computers of a molecule having a 
simpler coupling pattern, and eases the problem of un- 
wanted influences on spins. 

III. DESIGN OF TEST 

The proof of the pudding is in the eating: the 5- 
bit NMR computer to be described was tested on the 
Deutsch- Jozsa problem for functions of 4 bits Gj , in the 
form described in [[ll| , modified for efficiency with NMR 
as described by Jones and Mosca |Q. (A recent sim- 
plification |l3f| , unused here, would permit working with 
functions of 5 bits.) The problem is to decide whether a 
function program selected from a set of possible programs 
computes one kind of function or another. Specifically, 
the problem is to distinguish programs for balanced func- 
tions from programs for constant functions, where the 
functions are from {0, l} 4 to {0, 1}. (A function is con- 
stant if its value is independent of its argument, and is 
called balanced if the value for half the arguments is 1 
while the value is for the other half.) The test actu- 
ally made was to distinguish between programs for one 
constant and one balanced function, defined as follows: 

fo(x) d ^ (1) 

and 

f b (x) d = x x ®x 2 <Sx 3 ®x 4 (2) 

for all x, where x = (x\, X2, ^3,^4), and "©" is addition 
modulo 2. Also, several controlled-not (CNOT) gates 
were tested, along with a variety of 1-bit operators. The 
balanced function chosen, ft, has the nice property of 
being implementable also in classical reversible gates with 
no work bits. 

Used to solve this problem, a quantum computer is 
a resource used both to specify the function under test 



and to determine what it is. In order to separate these 
two uses, one can view the quantum computer as used 
alternately by a specifier of the function and a decision 
maker, two players of a game in which: (A) the decision 
maker prepares the starting state; (B) the specifier runs 
the function program;!!] and (C) the decision maker makes 
a measurement independent of the function program, and 
interprets the result to decide the function class. 

On an NMR quantum computer, (A) the decision 
maker starts a play by using r.f. pulses and magnetic-field 
gradients (independent of the function to be specified) to 
put the liquid sample in the pseudopure state having a 
density matrix with a traceless part proportional to 

Pl = f i6|ooooi)(ooooi| - -1 = m?i%i%ii4 - o 1 ( 3 ) 

in terms of the polarization operators 1^ = (^1 + Ikz) 
and if = (±1 - I kz ) usual to NMR |l0]||. Then the 
decision maker applies a unitary transform Ugo by use of 
a hard 90° y-pulse which for this particular state has the 
same effect as the Hadamard transform on each spin [^2) . 

Pi -* PO = U 90 PiU 90 

x (i 1+J -)(i 1 - J -)-i 1 - (4) 

(B) The specifier chooses a function / from one of the 
set of functions undergoing test, here /o or /{,, and runs 
the quantum version of a program to compute /; this 
program is a sequence of gates, each a unitary transfor- 
mation implemented by an r.f. pulse sequence. The total 
program implements a unitary transformation U(f), 

defined by its action on basis vectors \x, x§): 

U(f)\x, x 5 ) = \x, x 5 ®f(x)). (5) 

The transform U (/) produces the density matrix with 
traceless part proportional to pf. 



(C) The decision maker reads out the NMR spectrum 
which depends on pf. The spectrum differs according 
to whether / is balanced or constant, and thus tells the 
decision maker the function class, with only one func- 
tion evaluation, a large saving over classical computa- 
tion, which could require 9 evaluations for functions of 
four bits. 



1 These two moves must be iterated for a classical computer, 
but not in the quantum solution of the Deutsch- Jozsa prob- 
lem, giving the quantum computer a large advantage over the 
classical computer. 
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In theory, for the case f — fo, U(f) is specified to 
be U(fo), which by Eqs. (1) and (5) turns out to be 
the identity matrix, so one should have pf = po. The 
spectrometer detects only the terms of the righthand side 
of Eq. (|j) that are linear in I x , so for a spectrometer 
adjusted to give an upward peak for I x , the resulting 
spectrum is in theory i ) which has, from left to right, 
positive peaks for spins 1 to 4 and a negative peak for 
spin 5. 

For the balanced function, f = fb (Eq. @), U(fb) is 
defined by U(fb)\x, x$) = \x, xi®x 2 ®X3®X4®x 5 ) . A 
unitary operator that is simpler to implement, that has 
the same effect on the fifth (value) bit, and that allows 
the distinction between constant and balanced functions 
is U(fb) defined by 

U(f b )\x, x 5 ) = \xi, Xi © X 2 , Xi © X 2 © X 3 , Xi © x 2 

© x 3 © x i ,x 1 © x 2 © x 3 © x 4 © x 5 ), (7) 

which we implemented by sequential application of the 
gates (CNOT)i2, (CNOT) 23 , (CNOT) 34 , and (CNOT) 45 . 
The spectrum calculated for the density matrix pi, ob- 
tained by transforming po with U(fb) is (vide infra). 



IV. REALIZATION AND TEST OF A 5-BIT NMR 
COMPUTER 

A 5-bit NMR quantum computer requires a molecule 
having 5 spin-active nuclei, with long relaxation times. 
Large separation of resonance frequencies of the nuclei 
allows rapid selective control of the spins. For frequency 
separation, it is desirable to use different atomic species 
for different spins, which requires a multi-channel NMR 
spectrometer. Our NMR experiments were performed us- 
ing a BRUKER AVANCE 400 spectrometer with five in- 
dependent r.f. channels and a QXI probe (H,C-F,N). The 
lock coil was also used for deuterium decoupling utilizing 
a lock switch. A linear path of spin-spin couplings is suf- 
ficient for all computations || . Given the availability of 
a 5-channel spectrometer, we chose as the "hardware" of 
our NMR quantum computing experiments the molecule 
BOC-( 13 C2- 15 N- 2 D2 -glycine)-fluoride which contains an 
isolated coupling network consisting of five nuclei, each 
having spin 1/2: the amide 1 H, the 15 N, the aliphatic 
13 C Q , the carbonyl 13 C, and the 19 F nuclear spin (see 
Fig. 1). For simplicity, we will refer to these spins (and 
the corresponding bits) as 1, 2, 3, 4 and 5, respectively. 
All spins are heteronuclear, except for C a and C which 
however have a relatively large chemical shift difference. 
The five-spin system is well isolated from the protons of 
the BOC protecting group which are separated by more 
than four chemical bonds. In addition, the deuterium 
spins (D) which are attached to C" can be fully decou- 
pled from the spins of interest using standard heteronu- 
clear decoupling techniques |IJ,[l5|,|l6| . The 



substance was synthesized starting from commer- 
cially available 13 C and 15 N labeled glycine (see Ap- 
pendix A) and was dissolved in deuterated dimethyl- 
sulfoxide (DMSO-Dg). NMR experiments were per- 
formed at a magnetic field of about 9.4 Tesla and a sam- 
ple temperature of 27° C. The experimentally determined 
T2 relaxation times for spins 1-5 were 250 ms, 490 ms, 
450 ms, 590 ms, and 260 ms, respectively. Resonance fre- 
quencies Vk and scalar coupling constants Jki are summa- 
rized in Table I. Except for the J23 coupling constant of 
13.5 Hz, the spin chain is connected by one-bond coupling 
constants Jk{k+i} larger than 60 Hz. In the multiple ro- 
tating frame (see Ref. |u| and Appendix B) the preces- 
sion frequency of each individual spin is 0, which consid- 
erably simplifies implementation, because only coupling 
terms need to be considered (and manipulated). 

TABLE I. Resonance frequencies Uk, chemical shifts 5k, 
one-bond coupling constants Jk{k+i}, an <A non-zero two-bond 
coupling constants Jk{k+2} of the used five-spin system. No 
resolved three- or four-bond coupling constants were observed. 



V\ = 


400, 133,001.6 Hz 


(*i = 


7.51 ppm) 


l>2 = 


40, 547, 895.3 Hz 


(82 = 


75.54 ppm) 


V3 = 


100,616,858.0 Hz 


(63 = 


41.05 ppm) 


1/4 = 


100,629,089.1 Hz 


($4 = 


162.61 ppm) 


^5 = 


376,510,545.5 Hz 


(5 5 = 


31.92 ppm) 


Jl2 = 


= 94.1 Hz 


J23 = 


13.5 Hz 


J34 ~ 


= 65.2 Hz 


J45 = 


366.0 Hz 


Jl3 = 


= 2.7 Hz 


J35 = 


67.7 Hz 



TABLE II. List of initial Cartesian product operator terms 
of pi (Eqs. (3) and (^)) that give rise to detectable signals for 
at least one of the functions fo or ft in the implemented ver- 
sion |ll| of the Deutsch-Jozsa algorithm (see also Appendices 
B and C). The propagators U(f ) = 1 and U(f b ) (Eq. (7)) 
transform po (Eq. (4)) to po and p b , respectively (Eq. (6)). 
Only the underlined terms of po and p;, which contain sin- 
gle transverse spin operators correspond to single quantum 
coherences that are detectable in an NMR experiment. 
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PO 
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— Ibx 


— Ibx 



3 



The experimental implementation of the propagator 
U(fo) corresponding to fo is trivial because by Eq. (||) 
the propagator is the unit operator, implemented by do- 
ing nothing. In contrast, the construction of the pulse 
sequence to implement the series of CNOT-gates that 
define the unitary transformation U{fb) of Eq. (0) for 
the balanced function fb (Eq. (g)) requires attention. 
The goal is to create robust pulse sequence elements that 
minimize the effects of experimental imperfections. The 
pulse sequence elements shown in Fig. 2 A-D were de- 
signed specifically for the coupling topology of our 5-spin 
system to implement the unitary operators corresponding, 
to (CNOTjia, (CNOT) 23 , (CNOT) 34 , and (CNOT) 45 .@ 
During these CNOT gates that act on two directly cou- 
pled spins k and I, only the couplings Jki are active, while 
the effect of all other couplings in th e spin sy stem are re- 
Figure 3 



focused by cyclic pulse sequences 10 
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shows schematically the pulse sequence actually used for 
the propagator U(fb) for the balanced function /;,; this 
sequence benefited from applying simple rules for pulse 
cancellation (see Appendix B.2). 

The NMR implementation of the Deutsch-Jozsa algo- 
rithm starts with the preparation of the pseudopure state 
Pi of Eq. (0). The preparation of such a pseudopure 
state by a single pulse sequence requires a non-unitary 
transformation of the thermal equilibrium density oper- 
ator [f7|,[l9l . This can be achieved using spatial averaging 

In the basis formed by 
pi can be expressed as 



jj],H or temporal averaging (?| 
Cartesian product operators [1 
a linear combination of 31 terms that only consist of z 
spin operators: 



Pi ^ ^ &n Inz ^ ^ Sri ^Imzlnz ~t~ ^ ^ Sn ^-^Iz^mz^n 
n=l m<n l<m<n 

+ s nSI kz Ilz IGhzhzhzhzhz, 



k<l<m<n 



(8) 



where s n = — 1 if n = 5 and s n = 1 otherwise. It 
is straightforward to create each of these terms from 
the thermal equilibrium density operator, using standard 
building blocks of high-resolution NMR . In princi- 
ple, temporal averaging could be realized by repeating 
steps (A)-(C) of the game for all 31 terms in Eq. (@) and 



2 During each pulse sequence shown in Fig. 2 only the cou- 
pling Jfe; is active which is required in order to implement 
(CNOT)^;. The effects of other non-zero couplings (see Table 
I) are effectively eliminated, except for J13 in the sequence 
implementing (CNOT) 4 5 (Fig. 2 D). Although it would be 
straightforward to remove also this coupling, this would re- 
quire additional pulses which can be avoided because in our 
spin system the coupling J13 =2.7 has a negligible effect dur- 
ing the relatively short duration A45 — l/(2J4s) = 1.39 ms of 
this gate. 



by summing up the resulting spectra. However, because 
currently available NMR spectrometers require a distinct 
experiment to detect each spin species ( 1 H, 15 N, 13 C and 
19 F) (see Appendix B.3), a total of 124 NMR experi- 
ments would be required for each function / in order to 
include all terms in the temporal averaging. A detailed 
analysis shows that of the 31 terms that constitute the 
pseudopure state pi, only the five linear terms Ikz and 
the four bilinear terms 2IkzI{k+i}z are transformed into 
detectable operators by the propagator U90 (to create po) 
followed by the propagators U(fo) or U(fb), as the case 
may be (see Table II and Appendix C). As pointed out 
previously g], preparing just the 

linear terms Ikz suffices in some cases of the Deutsch- 
Jozsa problem to distinguish constant from balanced 
functions, because in these cases a balanced function 
gives a vanishing signal for at least one of the input spins. 
However, in the presence of experimental imperfections, 
it is desirable to identify a balanced function based on 
the sign reversal of the signal of at least one of the input 
spins, rather than by the lack of a signal. For the spe- 
cial case of the balanced function fb that was chosen for 
this demonstration experiment, this can be achieved by 
including also the bilinear terms 2IkzI{k+i}z as starting 
operators (see Table II). 

Samples described by these linear and bilinear terms of 
Pi were prepared (see Appendix B.4) to demonstrate ex- 
perimental control of the five-spin system and to execute 
cases of the Deutsch-Jozsa algorithm. For each function 
(/o and fb) the following three sets of experiments were 
performed (see experimental spectra in Fig. 4). Set 1 
(first row of curves from the bottom in Fig. 4): prepa- 
ration of the linear terms Ikz (with algebraic signs as 
specified in Eq. (8) and Table II), application of C/90 and 
U(f), and detection of spin k for k = 1, . . . , 5; set 2 
(second row in Fig. 4): preparation of the bilinear terms 
2IkzI{ k +i}z (with algebraic signs as specified in Eq. (8) 
and Table II), application of Ugo and U(f), and detection 
of spin k for k = 1, . . . , 4; and set 3 (third row in 

Fig. 4): preparation of the bilinear terms 21{ k -i}zlkz 
(with 

algebraic signs as specified in Eq. (8) and Table II), 
application of Ugo and U(f), and detection of spin k for 
k = 2, . . . , 5. 

The observed spectra shown in Fig. 4 correspond 
closely to the theoretical predictions (see Table II). For 
the constant function fo, only the experiments of set 1 
yield detectable signals. For the balanced function fb, 
the experiments of set 1 only yield a detectable signal 
for spin 5, whereas for spins 1-4 detectable signals are 
only obtained in the experiments of set 2. As expected, 
only spurious signals are detected for the experiments of 
set 3. The amplitude of these spurious signals is typi- 
cally on the order of 4% compared to the full signals. As 
expected, all the signals of spins 1-4 are positive for the 
constant function whereas the signal of spin 4 is inverted 
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by the propagator U(fb) corresponding to the balanced 
function. For fb the signal amplitudes reach only be- 
tween 55% and 70% of the amplitudes found for / . This 
signal loss can be attributed mainly to relaxation and 
experimental imperfections during the sequence that im- 
plements U(fb) (Fig. 3), which has an overall duration of 
51.4 ms. 

Through combined synthetic, analytic, and spectro- 
scopic work, a five-bit NMR quantum computer was built 
and shown to implement superposition, quantum inter- 
ference, and designed unitary transformations. Although 
obstacles had to be overcome, none were fundamental, 
and quantum computers with more than five bits will 
be built. Lots of interesting questions have been raised 
for future work pertaining to the constraints and oppor- 
tunities for linking molecular architecture, spectrometer 
design, and algorithms for NMR quantum computing. 
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APPENDIX A: SYNTHESIS OF MOLECULE 

We purchased 250 mg of 13 C2- 15 N-glycine from Martek 
Biosciences Corporation, 6480 Dobbin Road, Columbia, 
Maryland 21045. The labeled glycine was fully deuter- 
ated by treatment with a solution of NaOD in D2O 
at 140°C. The product was dissolved in water for re- 
protonation while retaining the deuterium atoms in 
alpha-position. The resulting 13 C2- 15 N- 2 D2 -glycine was 
protected in a standard reaction with di-ierf.-butyl- 
dicarbonate (BOC-anhydride) as reagent (O. Keller, 
W. E. Keller, G. van Look and G. Wersin, Org. Synth. 63, 
160 (1985)). Finally the carboxylic acid was converted 
by cyanuric fluoride into the desired acyl fluoride: BOC- 
( 13 C 2 - 15 N- 2 D^-glycine)-fluoride (L. A. Carpino, E. M. E. 
Mansour and D. Sadat- Aalaee, J. Org. Chem. 56, 2611 
(1991)). The substance dissolved in DMSO-D 6 at room 
temperature shows NMR spectra that weaken with a 
half-life of about a week, indicative of reactions not yet 
determined. The solution was stable during storage at a 
temperature of —30° C. 



APPENDIX B: NMR PULSE SEQUENCES 

For the preparation of the elements of a pseudopure 
state and the implementation of quantum gates, robust 
r.f. pulse sequences are desirable. Pulse-sequence param- 
eters with negligible experimental errors are the dura- 
tions of r.f. pulses and of delays. In addition, the phases 
of r.f. pulses and of the receiver can be controlled with 
negligible errors. The most important experimental im- 
perfections arc r.f. amplitude errors that result from mis- 
calibrations and from the r.f. field inhomogeneity created 
by the r.f. coils. In addition to the use of compensating 
schemes, such as super cycles and composite pulses [ jltfl , 
experimental imperfections can be reduced by designing 
pulse sequence elements with a minimum number of r.f. 
pulses. For example, pulses to refocus frequency offset 
terms in homonuclear spin systems with different chemi- 
cal shifts can be eliminated by implementing the exper- 
iments in the multiple-rotating frame in which the pre- 
cession frequency of each individual spin is (see section 
B.l). More generally, pulses can often be eliminated or 
replaced by phase adjustments with negligible errors (see 
section B.2). For the available spectrometer, the experi- 
mental pulse parameters are summarized in section B.3. 
The preparation of the elements of the pseudopure state 
Pi is discussed in section B.4. 

B.l. Implementation of experiments in the multiple 
rotating frame 

For heteronuclcar spins with resonance frequencies i>k 
and vi in the laboratory frame, the spins are irradiated 
on-resonance and the observed signals are demodulated 
by the determined resonance frequencies. If only a single 
r.f. channel is available for several homonuclear spins, on- 
resonance irradiation of several homonuclear spins can 
be achieved using phase-modulation of the r.f. pulses. 
The reference phase of each pulse applied to spin k must 
be adjusted such that it matches the desired phase in 
the corresponding rotating frame (vide infra). In addi- 
tion, the phases of the detected signals need to be cor- 
rected for the relative phases that have been acquired 
by the respective rotating frames during the course of 
the experiment. In our case with the two homonuclear 
spins C a (spin #3) and C (spin #4), the transmitter 
frequency of the carbon r.f. channel was set to the C a 
resonance frequency. In order to simplify the combina- 
tion of different quantum gates, the durations of the pulse 
sequences for each gate were chosen to be integer multi- 
ples of A = 1/|^3 — v±\ = 81.75 /is. Hence, the rotating 
frames are aligned at the end of each gate. 
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B.2. Simplifying pulse sequences 

Some quantum gates, such as (CNOT)^;, require z ro- 
tations of individual spins which can be implemented us- 
ing composite r.f. pulses [p2|. However, these pulses can 
be avoided if z rotations (by angle <p) are implemented 
by a corresponding negative rotation of the respective ro- 
tating frame of reference. In practice, this results in an 
additional phase shift (by angle — cp) of all following r.f. 
pulses that are applied to this spin and of the receiver 
phase of this spin. Furthermore, 180°, pulses (with arbi- 
trary phase #) are required in some cases to refocus the 
evolution due to J couplings. In order to undo the ro- 
tation caused by these pulses, additional 180°, or 180° ^ 
pulses are often needed at the beginning or at the end of 
these quantum gates. An appropriate choice of the posi- 
tion and phase d of these pulses often makes it possible 
to cancel two pulses from adjacent gates (e.g., 180° and 
ISO .;,,) or to absorb a 180° pulse into the phase of an ad- 
jacent 90° pulse (e.g., a 180° pulse preceded or followed 
by a QQ°_ X pulse is equivalent to a single 90° pulse). 

Even if r.f. pulses cannot be completely eliminated, the 
accumulation of small flip angle errors can be avoided 
by a proper choice of pulse phases which is common 
practice in the design of modern NMR multiple pulse 
sequence s [|10||. For example, the so-called MLEV-4 ex- 
pansion @|l|@ lSO^lSO^lSO^lSO . (used here, e.g., 
for spin 5 decoupling during spin 3- and spin 4-selective 
90° pulses) is preferable to lSO^lSO^lSO"^ ^ or to 
180°180°180°180°. 

B.3. Experimental pulse parameters 

Due to their large frequency separation, selective 
pulses for spins 1 ( 1 H), 2 ( 15 N) and 5 ( 19 F) could be 
implemented by simple square pulses. The durations of 
90° pulses were 8.85 [is, 41 /is and 11.75 [is, respectively. 
For spins 3 ( 13 C Q ) and 4 ( 13 C) the following shaped 
pulses with minimal durations and optimal selectivity 
were chosen based on numerical simulations and exper- 
imental optimizations: 90° pulses were implemented as 
e-SNOB pulses [||, not for the usual 270°, but for a 90° 

rotation with a duration of 224 [is; selective 180° pulses 
were implemented as Gaussian pulses [|24j with a duration 
of 250 [is and a truncation level of 20%. The application 
of these shaped e-SNOB and Gaussian pulses on C a has 
a nonresonant effect [^5| on C which corresponds to ex- 
perimentally determined z rotations of tp e = —4° and 
ip g = —18°, respectively. Conversely, a shaped e-SNOB 
pulse applied to C leads to a z rotation of — <p e for C a . 
In all experiments these phase shifts were taken into ac- 
count by adjusting the phases of the following pulse and 
the receiver phases (see Fig. 3). (Note that the phases of 
the two selective Gaussian 180° pulses applied to spin 3 
in the period T45 is not corrected because their absolute 



phases are arbitrary, c.f. Appendix B.2.) During the spin 
3- or spin 4-selective 180° pulses, the evolution due to the 
strong J35 and J45 couplings is automatically refocused. 
As this is not the case for spin 3- or spin 4-selective 90° 
pulses, spin 5 was actively decoupled during these pulses 
(see Fig. 3). 

As commercial high-resolution NMR spectrometers are 
commonly not equipped with multiple receivers, it was 
not possible to simultaneously detect the signals of differ- 
ent spin species. Moreover, the application of any given 
pulse sequence required four different pulse programs be- 
cause the routing of the r.f. channels (for the creation of 
1 H, 15 N, 13 C, 19 F and 2 D pulses) depends on the detected 
spin species ( 1 H, 15 N, 13 C or 19 F). Due to this techni- 
cal limitation, each spin species had to be detected in a 
separate experiment for every term of the initial density 
operator pi. However, this made it possible to use stan- 
dard heteronuclear decoupling techniques to simplify the 
detected signals and to significantly increase the signal- 
to-noise ratio of the experiments. 

During spin 1 detection, spins 2 and 3 were decoupled 
with an r.f. amplitude v r f = "/B r f /(2ir) of 0.6 kHz and 

0. 4 kHz, respectively. During spin 2 detection, spins 1 
and 3 were decoupled with an r.f. amplitude of 2.3 kHz 
and 0.4 kHz, respectively. During spin 3 detection, spins 

1, 2 and 5 were decoupled with an r.f. amplitude of 2.3 
kHz, 0.6 kHz and 2.0 kHz, respectively, and during spin 

4 detection, spin 5 was decoupled with an r.f. amplitude 
of 2.0 kHz. In all these cases, the WALTZ-16 decou- 
pling sequence 0Jl5|,[l6| was used. In principle, also 
the J34 coupling could be effectively eliminated during 
detection of spin 3 or spin 4 using time-shared decou- 
pling. However, this was not possible with our experi- 
mental setup because more than five separate r.f. chan- 
nels would have been required. From the resulting dou- 
blets (with splitting J34) an apparent singlet was created 
by merging the two doublet components |26|] . During spin 

5 detection, spins 3 and 4 were simultaneously decoupled 
using a double-selective G3-MLEV sequence |]l4] , ^5| , ^6| 
with an r.f. amplitude of 6 kHz. In addition, during all 
experiments deuterium decoupling was applied using a 
WALTZ-16 sequence with v r f — 0.5 kHz. In order to ap- 
proximate a constant sample temperature of about 27° C 
in spite of the additional sample heating effected by the 
decoupling sequences, 16 dummy scans were used prior 
to signal acquisition of spins 1 and 5, whereas 4 dummy 
scans were used prior to signal acquisition of spins 2, 3, 
and 4. Nevertheless, the linewidths of the experimental 
signals shown in Fig. 4 were slightly increased by resid- 
ual sample heating effects and imperfections of the de- 
coupling sequences. 



B.4. Pulse sequences for the preparation of the 
terms of p; and po 

In order to improve the sensitivity of the experiments 
and to filter out signals from impurities in the sample, in- 
dividual Cartesian product operator terms of pi were cre- 
ated using sequential INEPT transfer steps [21J starting 
from 1 H magnetization, corresponding to the operator 
l\ z . The term l\ z was prepared from the thermal equilib- 
rium density operator by applying spin 2, 3, 4, and 5 se- 
lective 90° pulses followed by a pulsed field gradient of the 
static magnetic field. An X-filter element 27 was used 
to select 1 H spins that are coupled to 15 N. For the prepa- 
ration of other terms of pi (starting from I\ z ), the phase 
4> a of the first 90° pulse applied to spin 1 was subject to 
a two-step phase cycle. In addition, the phase 4>b of the 
90° pulses for the implementation of C/ (see Eq. (4)) was 
also subject to an independent phase cycle. Overall, this 
resulted in a four-step phase cycle with the pulse phases 
a = {0°, 180°, 0°, 180°}, 4> b = {90°, 90°, 270°, 270°} and 
the relative receiver phases rcc = {0°, 180°, 0°, 180°}. 



APPENDIX C: DENSITY OPERATOR TERMS 

For all 31 Cartesian product operator terms in pi (Eq. 
(0)), the corresponding terms in po and pb are summa- 
rized in Table III. The transformation po U -Ji> p b of the 
unitary operator corresponding to the balanced function 
fb is composed of four consecutive unitary transforma- 
tions corresponding to (CNOT)^; quantum gates: 

(CNOT) 12 , (CNOT) 23 „ (CNOT) 34 ,„ (CNOT) 45 _ 
Po — > P — > P — ► P — > Pb- 

The transformations of the individual (CNOT)m gates 
can be derived using the following rules [0 : 



[5] 
[6] 
[<] 



'kx 


(CNOT) fc 


^^kx^-lx ) 




(CNOT) fc 


j 

1 kx ; 


h.r 


(CNOT) fc 





The terms of the intermediate operators p' , p" and p" 
are also given in Table III for completeness. 
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TABLE III. Terms of operators. 



Po 



Pb 



hz 
hz 
hz 

hz 

-hz 

2hzhz 

2hzhz 

Zhzhz 
—2hzhz 

2hzhz 

2hzhz 
—2hzhz 

%hzhz 
—Zhzhz 
—Zhzhz 

^hzhzhz 
^hzhzhz 
—^hzhzhz 
^hzhzhz 

-4hzhzhz 
-4hzhzhz 

^hzhzhz 
—^hzhzhz 
—^hzhzhz 
—^hzhzhz 

Shzhzhzhz 
Shzhzhzhz 
Shzhzhzhz 
—8h z h z hzhz 
Shzhzhzhz 

-lGh z hzhzhzhz 



hx 
hx 
hx 
hx 
~hx 

2hxhx 

2hxhx 

2hxhx 
—2hxhx 

2hxhx 

2>hxhx 
—^hxhx 

2hxhx 
—2hxhx 
—2hxhx 

4hxhxhx 

^hxhxhx 
—4I\xhxhx 

4hxhxhx 
4-hxhxhx 
—^hxhxhx 

4hxhxhx 
—4hxhxhx 
—^hxhxhx 
—^hxhxhx 

Shxhxhxhx 
Shxhxhxhx 
Shxhxhxhx 
Shxhxhxhx 
Shxhxhxhx 

-lGh x hxhxhxhx 



Zhxhx 
hx 
hx 
hx 
~hx 

hx 
4I\xhxhx 
4-hxhxhx 
—4I\xhxhx 
2hxhx 
2hxhx 
—2hxhx 
%hxhx 
—Zhxhx 
—^hxhx 

Zhxhx 
^hxhx 
^hxhx 
Shxhxhxhx 
Shxhxhxhx 
Shxhxhxhx 
4-hxhxhx 
—4hxhxhx 
—^hxhxhx 
—^hxhxhx 

^hxhxhx 
—4I\xhxhx 
—4Ii x hxhx 
-I6h x hxhxhxhx 
Shxhxhxhx 

Shxhxhxhx 



4hxhxhx 
2hxhx 
hx 
hx 
~hx 

hx 
^hxhx 
Shxhxhxhx 
Shxhxhxhx 
hx 
^hxhxhx 
—^hxhxhx 
%hxhx 
—Zhxhx 
—2hxhx 

Zhxhx 
^hxhx 
—Ihxhx 
4hxhxhx 
—4I\xhxhx 
-lfthxhxhxhxhx 
Zhxhx 
^hxhx 
Shxhxhxhx 
—^hxhxhx 

^hxhxhx 
—4I\xhxhx 
^hxhxhx 
Shxhxhxhx 
—^hxhxhx 

Shxhxhxhx 



Shxhxhxhx 
^hxhxhx 
2hxhx 
hx 
~hx 

hx 
2>hxhx 
4I\ x hxhx 
-Whxhxhxhxhx 
hx 
2hxhx 
Shxhxhxhx 
hx 

—4hxhxhx 
—^hxhx 

4-hxhxhx 
2>hxhx 
^hxhx 
4-hxhxhx 
—4I\xhxhx 
Shxhxhxhx 
Zhxhx 
—%hxhx 
—^hxhxhx 
—2hxhx 

2hxhx 
Shxhxhxhx 

^hxhxhx 
Shxhxhxhx 
—^hxhxhx 

—4I\xhxhx 



^hxhxhxhxhx 
Shxhxhxhx 
^hxhxhx 
^hxhx 
~hx 

hx 
2hxhx 
4I\ x hxhx 
Shxhxhxhx 
hx 
2hxhx 
—^hxhxhx 
hx 
—Zhxhx 
—hx 

Shxhxhxhx 
^hxhxhx 
^hxhx 
Shxhxhxhx 
—4I\xhxhx 
—Shxhxhxhx 
4-hxhxhx 
—Zhxhx 
—^hxhxhx 
—2>hxhx 

2hxhx 
—4hxhxhx 

—Zhxhx 
—^hxhxhx 

—%hxhx 

—4I\xhxhx 



Note: detectable terms of po and pt are underlined. 



